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(54) Communications system having pilot signals transmitted over frequency divided channels 



(57) In a communications system signals, which include predetermined pilot signals, are transmitted over 
selected channels of a number of frequency-divided channels, such as in an M-QAM system. The receiver 
comprises means for receiving predetermined pilot signals simultaneously on selected (e.g. outer) channels of 
the plurality of channels and for calculating (at 33) channel gain values for said channels, first computation 
means for calculating frequency domain characteristics (r.m.s. delays) of the channels representative of the 
relative behaviour of the channels, based on the reception of pilot signals received simultaneously on said 
selected channels, second computation means for interpolating or extrapolating in the frequency domain the 
behaviour of channels on which no pilot signal is received in a given timeslot based on the behaviour of 
channels on which pilot signals are received in that timeslot to provide channel gain values for the remaining 
channels and third computation means responsive to said channel gain values, for demodulating the signals 
received on the channels on which no pilot signals are received. 
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Communications System Having Pilot Signals Transmitted Over Frequency 
Divided Channels 

Field of the Invention 

5 

This invention relates to a communications systems where signals, 
which include predetermined pilot signals, are transmitted over selected 
channels of a number of frequency-divided channels, such as in an M-QAM 
system. 

10 

Background of the Invention 

One of the problems associated with highly spectrum efficient digital 
transmission over land mobile channels is the frequency selectivity of the 

15 channel together with its Rayleigh fading characteristics. To overcome this 
problem, two solutions have been used: adaptive equalization or frequency 
division multiplexing (FDM). This invention is primarily but not exclusively 
related to the latter. With this technique, the entire channel (25 KHz for 
example) is subdivided into subchannels. The subchannels are narrow 

20 enough so that the transmission media over this reduced bandwidth can be 
considered as non-frequency selective or having frequency flat fading (F 3 ) 
characteristics. This technique is denoted Multi-Quadrature Amplitude 
Modulation or M-QAM. 

In every subchannel, a sequence of predetermined pilots is embedded 

25 in the information in order to compensate for the flat fading introduced by 
the channel. 

A block diagram of a baseband implementation of such a system is 
depicted in Fig. 1 for the particular case of 4-QAM. The input bit stream at a 
rate of R bits/sec is subdivided into 4 bit streams at rates R/4 bits/sec. A QAM 

30 modulator transforms every 4 bits in the subchannel into a single QAM 

symbol drawn from a 16QAM constellation (see for example, Proakis, Digital 
Communications, Mc-Graw Hill). Known symbols termed pilots are inserted 
in every one of the subchannels at predetermined times NT where T is the 
symbol rate. The resulting symbol rate for each subchannel in this particular 

35 scheme is equal to 4 Ksymbols/sec. The resulting QAM signal is then filtered 
with a pulse shaping filter (PSF). Each subchannel signal AM modulates an 
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oscillator centered at the subchannel frequency offset (offset w.r.t. the center 
frequency of the channel). The resulting signals are then added and translated 
to the actual transmission frequency. The purpose of the PSF in the 
subchannel branches is twofold. It limits the spectrum of the subchannel 
signals in order to avoid interchannel interference and the emission (splatter) 
into adjacent channels. The spectrum at the output of the transmitter is 
schematically depicted in Fig. 2. 

The signal is transmitted through the land mobile channel. This 
channel can be characterized as a multiple ray channel where every ray is 
independently Rayleigh faded having a gain equal to cr,-. A channel model is 
depicted in Fig. 3. The channel is characterized by its rms delay spread defined 
as 



(This notation assumes that the average delay is equal to zero or has been 
subtracted from the individual delays). 

Normally, when the rms delay of the channel is not too large with 
respect to the inverse of the bandwidth of the channel the multiple ray 
model can be simplified to a two ray model having the same rms delay 
spread. This simplified channel is depicted in Fig. 4 where t 0 is the average 
channel delay. Each channel is individually Rayleigh faded. The spectrum of 
the Rayleigh fading signal is characterized by its maximum Doppler frequency 
f d given as 

fd =1 c fc 

where v is the vehicle speed, c is the speed of light and f c is the actual 
transmission frequency. For 1.5 GHz and at a speed of some 100 km/h, f d is 
equal to 140 Hz. This Rayleigh fading can be modeled in baseband as a 
complex Gaussian process with in-phase and quadrature components having 
power spectral densities equal to 

^(/)=[i-(///,) 2 r i/2 

for \f\<f d and zero otherwise. This process is a relatively slow varying process 
with respect to the symbol rate (140 Hz vs. 4 KHz). 

The signal at the receiver is contaminated by white Gaussian noise. 
With an oscillator centered at each subchannel frequency, the signal is split 
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into its four subchannel components. A low pass filter in each subchannel 
rejects the signal corresponding to unwanted subchannels as well as noise. 
This signal is suitably sampled at the symbol rate. The amplitude and phase of 
the signal in each subchannel needs to be corrected for the variations 
5 introduced by the slow varying Rayleigh fading process. This is done in the 
following way: the gain at the pilot instants is determined by dividing the 
received signal at the pilot instant by the known pilot symbol. Once gains are 
known at the pilot instants and since the Rayleigh fading is slow compared to 
the symbol rate, interpolation can be done in this time domain to derive the 

10 gains at every symbol instant. The samples at the symbol instant are then 
divided by the estimated gain of the channel through this interpolation 
process and compared to a set of predetermined thresholds to determine the 
actual QAM symbol transmitted. The symbol is then converted to bits. The 
bits of the four subchannels are then combined to produce the original 

15 transmitted bit stream. 

Since the Rayleigh fading process has a certain limited bandwidth, 
according to the Nyquist criteria, in order to be able to recover this process, 
samples of this process need to be measured at at least twice the maximum 
frequency of the process. In the example above, f e was equal to 140 Hz and 

20 pilots ideally need to be included in each stream at a rate of at least 280 Hz, or 
280 pilots /sec. This assumes perfect frequency synchronization between 
transmitter and receiver. Although an automatic frequency control (AFC) 
circuit is normally implemented in such a system, some residual frequency 
offset between transmitter and receiver is normally present after AFC 

25 correction. In this case, the low pass Rayleigh fading process as seen by the 
receiver will not be centered around zero but around this frequency offset, 
increasing the perceived bandwidth at the receiver. To compensate for that, 
more pilots need to be introduced. Let us assume a residual offset of some 200 
Hz. The maximum frequency of the combined process (frequency offset plus 

30 Rayleigh fading) will be equal to 340 Hz. Again, using the Nyquist criteria, 
pilots need to be introduced at a rate larger that 680 Hz. A more real situation 
where interpolation cannot be accomplished with an ideal brick wall filter, 
some 800 pilots /sec will be more desirable. This means that one out of every 5 
symbols (800 pilots/sec out of 4 ksymbol/sec) in each subchannel is 'wasted' 

35 in this correction process, that is 20% of the information transmitted through 
the channel. 
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A possible prior art way to reduce this overhead is the following. Let 
the four subchannels be denoted as left/right (L/R) inner/outer (I/O), i.e., LO, 
LI, RI, RO. Pilots are introduced alternatively in the inner and outer channel 
in a way that will be explained below. The information in the channel is send 
5 in time slots using 4-16QAM modulation as described above. At the 
beginning of each slot, a sync pattern is used to aid in the time recovery 
process and AFC. The transmission in the outbound link is continuous, 
meaning that each slot is immediately followed by the next slot . The sync 
pattern consists of 3 predetermined symbols for each subchannel for a total of 

10 16 symbols. Now, instead of sending pilots in each subchannel as described 
above, they are distributed between inner and outer subchannels. The 
distribution of pilots in a slot is presented in Fig. 5. The subchannel is 
sufficiently narrow so mat one can assume flat fading (non-frequency 
selective) for each subchannel. In the time domain interpolation process, the 

15 gain at the pilot instants is determined for the inner subchannels and the 
outer subchannels. To find the gain for any other symbol time say in one 
outer channel, LO for example, the gain at the pilot instant in the inner 
channel (LI in this example) together with the gain in the outer channel are 
taken into consideration. This off-channel gain however cannot be used as if 

20 it were an on-channel gain since over the bandwidth of two subchannels, the 
fading cannot be assumed to be non-frequency selective. To alleviate the 
problem, in the time domain interpolation process a different weight is given 
to the off-channel gain at the pilot instants. This reduced weight was 
optimized through simulations and is constant independent of the rms delay 

25 present in the channel. This approach is not in anyway optimum. 

Summary of the Invention 

According to the present invention, a communications system is 
30 provided comprising a transmitter and a receiver, the transmitter comprising 
means for transmitting signals over a plurality of time slots on a plurality of 
frequency-divided channels (which in the preferred embodiment are sub- 
channels of a broader channel), said signals including predetermined pilot 
signals, wherein the pilot signals are transmitted simultaneously on selected 
35 channels of the plurality of channels and the receiver comprising: means for 
receiving predetermined pilot signals simultaneously on selected channels of 
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the plurality of channels and for calculating channel gain values for said 
channels, first computation means for calculating frequency domain 
characteristics of the channels representative of the relative behaviour of the 
channels, based on the reception of pilot signals received simultaneously on 
5 said selected channels, second computation means for interpolating or 

extrapolating in the frequency domain the behaviour of channels on which 
no pilot signal is received in a given timeslot based on the behaviour of 
channels on which pilot signals are received in that timeslot to provide 
channel gain values for the remaining channels and third computation 

10 means responsive to said channel gain values, for demodulating the signals 
received on the channels on which no pilot signals are received. 

In this manner, frequency domain interpolation or extrapolation is 
used to improve the receiver operation. The frequency domain 
characteristics are preferably the r.m.s. delays of the channels (i.e. the sub- 

15 channels). 

The invention can be employed in a system having pilot signals on all 
the channels, but it is particularly preferred that the transmitter transmits the 
pilot signals on those channels that are at the maximum separation in 
frequency. I.e., in the case of a number of sub-channels of a broader channel, 
20 the pilot signals are transmitted only (or substantially only) in the outermost 
sub-channels. This allows for interpolation between the frequencies (of the 
sub-channels) rather than extrapolation and provides better results. 

A preferred embodiment of the invention is now described, by way of 
example only, with reference to the drawings. 

25 

Brief Description of the Drawings 

Fig. 1 shows a block diagram of a typical prior art 16 QAM transmitter; 
Fig. 2 shows a schematic representation of the spectrum at the output 
30 of the transmitter of Fig. V, 

Fig. 3 shows a channel model for explanation of a typical radio channel; 
Fig.4 shows a simplified model of a radio channel; 
Fig. 5 shows a prior art time division multiplexed communication 
channel slot configuration; 
35 Fig. 6 shows a typical two-ray Rayleigh fading model for explanation of 

the invention; 
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Fig. 7 shows a block diagram of a receiver for explanation of the present 
invention; 

Fig. 8 shows the results of simulations carried out at different noise 

levels; 

Fig. 9 shows the proposed rms delay estimator as a function of the rms 
delay in the channel (or subchannel); 

Fig. 10 shows a new slot configuration in accordance with a preferred 
embodiment of the present invention; and 

Figs. 11, 12, 13 and 14 show different anticipated bit error rates for 
different signal-to-noise ratios in different environments. 

Detailed Description of a Preferred Embodiment 

In at least its preferred embodiment, the present invention requires 
that the RMS delay in the channel is known or is calculated. A preferred 
method of estimating the RMS delay is provided by the following 
mathematical derivation. 

Looking first at the typical two ray Rayleigh fading model shown in Fig. 
6 where the average delay has been ignored for simplicity. We assume that 
the first ray is delayed by an amount while the second one is delayed by - 
where is the rms delay in the channel. Both rays are then independently 
Rayleigh faded, x^t) and x 2 (t) represent the two complex Gaussian processes 
associated with the Rayleigh fading. 

The impulse response of the system can be expressed as 
g(t)=x 1 8(t-x)+x 2 5(t+T) 
Note that x, and x 2 are time dependent. 

Taking now the Fourier transform, a time varying frequency response 
is obtained thus: 

G, (/) = *i exp(-or) + x 2 exp(coT) 
A sync pattern is sent at the beginning of each slot. During the sync 
symbols, we can estimate the gain of the channel at four different frequencies, 
namely the subcarrier frequencies. Note that the time for all the symbols is 
identical. If we assume negligible noise in this estimation process, we may 
write 
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G(/ 0 ) = ^exp(-6) 0 T)+j: 2 exp(fi) 0 T) 

G(f 1 ) = x 1 cxp(-co l T)+x 2 exp(.(D 1 T) 

G(/ 2 ) = *i exp(-%T)+x 2 exp(o 2 T) 

G(/ 3 ) = x x exp(- o 3 t) + x 2 exp( a> 3 t) 
where f 0 ,f v f 2 ,f 3 af e the four subcarrier frequencies and <a, =2^ ; .In this 
particular case, these subcarrier frequencies are symmetrically spaced around 
zero and are related as follows. 

/o="3/ 2 
/i=-/ 2 
/ 3 = 3/ 2 

Substitution of this relation in the above equation yields: 

G(/„) = *, exp(3ffl 2 T)+ x 2 exp(-3w 2 T) 

G(/,) = x x exp(o 2 T)+ x 2 exp(-ffl 2 T) 

G(/ 2 ) = ^ exp(-£0 2 T) + x 2 exp(o 2 T) 

G(/ 3 ) = exp(-3o) 2 r)+x 2 exp(-3© 2 T) 
From this last set of equations, there exist several ways of determining the 
rms delay in the channel. All of them in some way rely on the fact that if the 
channel is frequency selective, the complex vectors representing the gain will 
not be collinear. A simple form is derived next. 
Let us define 

«02 = G(/ 0 )G*(/ 2 ) 
«i3 = G(/ 1 )G*(/ 3 ) 
aoo=G(/o)G'(/o) 
a^G^GV,) 

«22=G(/ 2 )G*(/ 2 ) 

a 33 = G(/ 3 )G*(/ 3 ) 

Next define 

A = a n + a 22 -a 02 -a, 3 
By direct substitution, it is easily shown that 

A = 2(|x 1 | 2 +|x 2 | 2 )(1-cos4o 2 t) 

If we denote now 

where E|jc,-| 2 is the signal energy in ray i, which in our model is identical for 
the two channels, then 

E(/l) = 4of(l-cos4(» 2 T) 
Now, if a) 2 T«lwe can approximate cosCx)=l-x 2 /2 and 
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E(A) = 32(j>fT 2 

that is E(A) is proportional to the square of the rms delay spread. 
Since a r is not really known, let us define 

/i = a 00 + a u + a 22 + a 33 
5 It can be easily shown that since E^jt^O (independent fade assumption) 
then 

We define the ratio 

E C*) A 2.2 

10 This ratio is independent of the signal energy in the rays and proportional to 
the square of the rms delay of the channel. 

The above derivation can be generalized to multiple rays as follows. 
Assume that we have n independently faded rays . The transfer function of 
the channel can be written as 

15 G(/)=5> t .expO-£0T ; ) 

We may assume that the clock recovery circuit in the system takes care of the 
average delay so that 

We define, as above, the set of variables a. We have that 

» 2 

a 02 = ex P(-y 4 %T ; ) + cross terms 

i=l 

n 2 

20 «i3 = XW expO'4%T,)+ cross terms 

«=i 

n 2 

a « = XW + cr oss terms 

i=l 

were the term 'cross terms' designates products of the form with 
We define, as above, 

^ = «11 + «22 ~ a 0 2 - a i3 

and 

25 /z = a 00 + a u + a 22 + a 33 

Since we assume independent fades for the rays, the E(x I oc*)=Ofor 
With this condition, it is easily shown that 

EU)=2^ of (l-a^a^.) 
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and 



where 



5 




2 V" of 



Assuming now that w 2 t, «1 for all i, we can rewrite 



10 




It has been shown that this quite simple procedure can help determine 
the rms delay in the channel. As mentioned, the procedure is based on the 
fact that the gains for each subchannel are not collinear. Other procedures can 

15 be devised that use this same fact. 

A simulation program was written to find out for a two ray model the 
validity of our derivation. A block diagram of the receiver is depicted in Fig. 
7. The receiver comprises a receive signal input 20 for receiving a digitized 
received signal from an R.F. stage (not shown). The received signal is 

20 simultaneously processed in four paths by being fed (in software) to four 
mixers 21-25 for injecting four subcarriers and in effect reducing the four 
components to zero i.f. One pulse shaping filter (PSF) 26 - 29 is provided in 
each path and the outputs of the PSFs are sampled by sampler 29 every T 
seconds (i.e. 1/T is the symbol rate). The samples are fed to storage and 

25 processing means 33 which stores the samples, calculates the channel gain 
values as described above, multiplies the samples by the gain values and 
performs threshold comparisons to provide the received bits, which are 
output as a received bit stream 34. 



30 mobile receiver can use in its delay estimation process the gains 

corresponding to the sync pattern in its own slot and the sync pattern 



Since, in the outbound channel a slot is followed by another slot, the 
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corresponding to the next slot, for a total of six sync symbols, and are 
computed for every sync symbol and their quotient is averaged over the six 
sync symbols. There is a possibility of averaging over more that one received 
slot to obtain more accurate results. 

In the above derivation, it is assumed that the noise associated with the 
gain estimation process is negligible. To see the effect of the noise, 
simulations for different cases have been performed. The estimation is 
averaged over six received slots. Fig. 8 shows the results obtained for as a 
function of E D /N 0 , the average energy per transmitted bit over the energy per 
unit of bandwidth of the noise process. The 'one sigma' points are also 
shown in the graph. It is seen that for Et,/N 0 exceeding some 20 dB, the noise 
has no effect on the estimator. 

Fig 9 presents values for as a function of the rms delay. The parabolic 
nature of the behavior is clearly seen. The small variance of is also evident. 

Based on the knowledge of the rms delay spread in the channel, it will 
be shown that by sending pilots in two of the subchannels, the gain in the 
other subchannels can be estimated through a straight forward frequency 
domain interpolation technique. One possible slot structure that uses pilots in 
two subchannels at a time is depicted in Fig. 10. The same technique can be 
applied to any other slot configuration that sends two pilots at a time, not 
only in the case of 4-QAM considered here but for a general M-QAM with 
M>2. 

Observe that, in this example, the pilots have been placed in the outer 
subchannels. They could have been placed in the inner subchannels or one 
inner and one outer, etc. This configuration however that uses 
'interpolation' in the frequency domain rather that 'extrapolation' offers 
better results. 

Interpolation is done in two steps: frequency domain and time domain. 
Since the pilots are known to the receiver, the gain at the pilot instants can be 
determined by dividing the received sample by the corresponding pilot 
symbol. In this case, we will not assume negligible noise for this gain 
estimation. 

Assume that the gains at the pilot instants for subchannels 0 and 3 can 
be expressed as 

§0=80 + 710 
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where & is the actual gain in subchannel i while 77, is the 'measurement 
noise'. The problem can be postulated as follows: Find a 0 and such that 

and 

r = E|ft-&| 2 -»min 

Direct substitution yields 

81-81= Oo(£o + Vo)+ 03(83 + %)-8i 
From the two ray model, we can write 

g 0 = exp(-j3%T)+Xj expO"3© 2 T) 
ft = x x exp(-j'o 2 T)+x 2 exp(;'ft) 2 T) 
g 2 = jc, exp(y3fi» 2 T)+ x 2 exp(-y3co 2 T) 
With the aid of this last equation, we can compute 
r = EH, -ft | 2 = 4(2e? + v 2 )+ <£(2o? + v 2 ) 
+4^0^ cos(6fi) 2 T)-4a 0 Cv cos(2cj 2 t) 
-403a 2 cos(4 co 2 T) + 2 a 2 

where 

v 2 = e|t/ 1 .| 2 fori = 1,2 

is the power of the noise associated with the gain measurement process. 

In order to find the minimum with respect to a 0 and , we 
differentiate and equate to zero 

— — = 0 => a 0 ( 2o v + v 2 ) +2a 3<^ cos(6© 2 t) = cos(2co 2 t) 

= 0 => 2ooO^ cos(6tu 2 T) + «3(2a 2 + v 2 ) = a 2 cos(4o 2 t) 

A system of two equations in two unknowns has to be solved. The solution to 
this system of equation is 

_ COS2ft> 2 T-COS6(P 2 TCOS4(P 2 T + ^" 1 COS2(P 2 T 

a ° sin 2 6eo 2 T + 2%~ l + x~ 2 

_ COS4(P 2 T-COS6fl; 2 TCOS2Q? 2 T+^~ 1 cos4o? 2 t 
^ ~ sin 2 6© 2 t+2^ -1 + x~ 2 

where we have defined 

lot 
X = -f 

which is related to the signal-to-noise ratio in estimating g 0 or g 3 . Optimal 
frequency domain interpolation coefficients have been derived for this case. 
From symmetry considerations, we can write the following equations for the 
two cases 
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It has been shown that from the sync part of the slot, the rms delay in 
the channel can be estimated. In the two outer channels, gains are estimated 
using the known pilots and sync symbols. Known sync symbols are also used 
5 in the inner channels to estimate the gain at the two ends of the slot (in the 
outbound case), while in the middle of the slot, frequency domain 
interpolation can be used. 

Once the gains are known at all the corresponding pilot instants, time 
domain interpolation can be done. To derive time domain interpolation 
10 coefficients, several techniques can be used. The simplest is to use a low pass 
interpolating filter with predetermined coefficients. 

Results for several typical environments characterized by different rms 
delay spreads obtained with this new technique are presented in Fig. 10 to 14 
and compared with the suboptimum technique described in the previous 
15 section. The center frequency of the system for all cases is assumed to be equal 
to 900 MHz in both cases. It can be observed there that at a 1\% Bit Error Rate, 
typical of a mobile communication system, an improvement of at least 1 dB 
can be obtained using this new technique compared to the old one. 

One possible variant of this new technique would be to assume a 
20 certain constant rms delay in the channel. In this case, frequency domain and 
time domain interpolation equations can be combined into a single set of 
equation for the inner channels. 

Although the above examples have emphasized the 4-16QAM 
technique, it should be dear the this same technique can be used for other 
25 multiple modulation schemes using this type of FDM technique with pilot 
insertion to correct for the unknown gain introduced by the frequency 
selective Rayleigh fading channel. As examples, all classes of M-NQAM or M- 
NPSK van be considered where M stands for the number of FDM 
subchannels and N for the number of points in the quadrature amplitude 
30 modulation (QAM) or phase shift keying (PSK). 
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Claims * 

1. A communications system comprising a transmitter and a receiver, 

5 the transmitter comprising means for transmitting signals over a plurality of 
time slots on a plurality of frequency-divided channels, said signals including 
predetermined pilot signals, wherein the pilot signals are transmitted 
simultaneously on selected channels of the plurality of channels and 

10 the receiver comprising: 

means for receiving predetermined pilot signals simultaneously on 
selected channels of the plurality of channels and for calculating 
channel gain values for said channels, 

15 

first computation means for calculating frequency domain 
characteristics of the channels representative of the relative behaviour 
of the channels, based on the reception of pilot signals received 
simultaneously on said selected channels, 

20 

second computation means for interpolating or extrapolating in the 
frequency domain the behaviour of channels on which no pilot signal 
is received in a given timeslot based on the behaviour of channels on 
which pilot signals are received in that timeslot to provide channel 
25 gain values for the remaining channels and 

third computation means responsive to said channel gain values, for 
demodulating the signals received on the channels on which no pilot 
signals are received. 

30 

2. A communications system according to claim 1, wherein the 
transmitter transmits the pilot signals on selected channels only and no pilot 
signals on the other channels. 
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3. A communications system according to claim 2, wherein, of the 
plurality of channels, the selected channels are at the maximum separation in 
frequency. 

5 4. A communications system according to claim 1, 2 or 3, wherein the 
frequency domain characteristics are the r.m.s. delays of the frequency divided 
channels. 

5. A communications receiver for receiving signals over a plurality of 
10 time slots on a plurality of frequency-divided channels, comprising: 

means for receiving predetermined pilot signals simultaneously on 
selected channels of the plurality of channels and for calculating 
channel gain values for said channels, 

15 

first computation means for calculating a frequency domain 
characteristic of the channels representative of the relative behaviour 
of the channels, based on the reception of pilot signals received 
simultaneously on said selected channels, 

20 

second computation means for interpolating or extrapolating in the 
frequency domain the behaviour of channels on which no pilot signal 
is received in a given timeslot based on the behaviour of channels on 
which pilot signals are received in that timeslot to provide channel 
25 gain values for the remaining channels and 

third computation means responsive to said channel gain values, for 
demodulating the signals received on the channels on which no pilot 
signals are received. 

30 

6. A communications receiver according to claim 5 further comprising 
fourth computation means for interpolating in the time domain the 
behaviour of all the channels based on the behaviour of channels on which 
sequential pilot signals are received at sequential times and based on 

35 frequency domain characteristics of the channels at said sequential times of 
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reception of the sequential pilot signals, thereby calculating channel gain 
values for all the channels in all time slots. 
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